      Subroutine Calc_Pressure(tRho, tP)
      include "param.fi"
      include "common.fi"
      REAL tRho, tP
      tP = CSON*CSON*RDENSITY/cGAMMA*((tRho/RDENSITY)**cGAMMA-1.0)
      End Subroutine

C==================================================================
      Subroutine PressureField
      include "param.fi"
      include "common.fi"

      Integer I
      Do I = 1, NTOT3T
        Call Calc_Pressure(Rho(I), P(I))
      Enddo

      End Subroutine
